Log-normal Distribution of Level Curvatures in the Localized Regime: analytical 

verification 
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We study numerically and analytically the moments of the 
dimensionless level curvature for one-dimensional disordered 
rings of the circumference L pierced by a magnetic flux (f>. 
The negative moments of the curvature distribution can be 
evaluated analytically in the extreme localization limit. The 
ensuing small curvature asymptotics of the corresponding dis- 
tribution has a "log- normal" behavior. Numerically studied 
positive moments show differences from other log-normally 
distributed quantities. 



Recently there was a considerable interest in various 
statistical characteristics of spectra of disordered and 
chaotic quantum systems, see e.g. Q. One particularly 
interesting issue is the so-called "level response statis- 
tics" characterizing sensitivity of individual energy levels 
E n (X) with respect to perturbation in an external param- 
eter A. The role of such a parameter can be played by 
e.g. an external electric or magnetic field, the strength 
and form of the potential, or any other appropriate pa- 
rameter of different nature on which the system Hamil- 
tonian is dependent. A convenient quantitative measure 
of level sensitivity is provided by a set of first and sec- 
ond derivatives dE n {\) / d\ and d 2 E n (X)/d\ 2 , known as 
"level velocities" and "level curvatures", respectively. 

The statistics of these quantities are mostly studied 
for the systems with completely "ergodic" eigenfunctions 
covering randomly, but uniformly all the available phase 
space and showing no specific internal structure. For sys- 
tems of this type most of their statistical characteristics 
are known to be universal, i.e. independent on particular 
microscopic details, and adequately described by ensem- 
bles of large random matrices of particular global sym- 
metry [0J^]. The same universality class comprises also 
weakly disordered (metallic) systems as long as effects of 
Anderson localization are negligible |l| [| . 

The form of the level curvature distribution typical for 
random matrices of various symmetry classes was guessed 
by Zakrzewski and Delande ||] on the basis of the numer- 
ical data. The analytical derivation of the Zakrzewski- 
Delande distribution is due to von Oppen || and Fyo- 
dorov and Sommers || . 

The effects of eigenfunction localization are expected 
to modify the level curvature statistics drastically. In 
contrast to the "random matrix" regime it is obvious 
that the level response should be less universal depending 
on the nature of perturbation. Such a non-universality 



manifests itself already in the form of perturbative cor- 
rections to the curvature distribution due to localization 
effects §. 

Actually, the sensitivity of energy spectrum to a 
change in boundary conditions was suggested long ago 
by Thouless as a measure of system conductance ||. To 
quantify this statement let us consider a one-dimensional 
sample closed to form a ring of circumference 2L encir- 
cling the Aharonov-Bohm flux cj) (measured in units of 
the flux quanta (f>Q = ch/e). The wavefunction in such 
geometry acquires the phase when going around the flux: 



m(x + L) =e i H(x-L). 



(1) 



The curvature is given by the second-order perturbation 
theory: 



K n = 



d 2 E n {4>) 



(m\P x \ n) 



(2) 



up to a constant shift. Here | m) and E m (0) are eigen- 
vectors and eigenvalues of the Hamiltonian at zero flux 
and P x is the momentum operator. Using a similarity 
of this expression to the conductance given by the Kubo 
formula, Thouless argued that the "typical" dimension- 
less level curvature K/A (with K measured, e.g. by the 
widths of the curvature distribution) is proportional to 
the dimcnsioneless conductance of the sample. 

Thouless original qualitative arguments seemed to be 
controversial and also did not take into account strong 
correlations between energy levels of a disordered sys- 
tem. The curvature-conductance relation gave rise to a 
lot of discussion and even claimed to be incorrect ||. 
The problem was reconsidered in much details recently 
ptj| , ^2 11 1, the results favouring validity of the Thouless 
idea in the metallic regime. Actually, for a disordered 
system in good metallic regime the Thouless relation, 
including the universal proportionality coefficient, was 
derived explicitly [|l3j in the course of analytical verifi- 
cation of the Zakrzewski-Delande distribution. However, 
the perturbative localization corrections to both quanti- 
ties are different Q. 

As to the issue of the level curvatures in systems with 
strongly localized eigenstates, our present understanding 
is based mainly on the results of numerical simulations 
p^| , p^| -p^]. It turns out that the curvature distribution 
is close to log-normal one for strong enough localization 
and has quite a non-trivial form in the vicinity of the 
Anderson transition |l5). Qualitative origin of the log- 
normal distribution can be inferred from Eq.(^) as sug- 
gested in pa]. Indeed, it is natural to assume that the 
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absolute value of the curvature \K\ is proportional to the 
product of amplitudes of a typical wavefunction at op- 
posite edges ±L of the sample: \K\ oc |^(— L)^{L)\ ~ 
cxp — (constL/£). The log- normal distribution follows 
if one assumes that the inverse localization length 
shows Gaussian eigenstate-to-eigenstate fluctuations. 

Whatever attractive and transparent is this simple ar- 
gumentation, a close inspection shows that it should be 
taken with caution. To this end let us recall that the 
distribution of the quantity v = \^(—L)^(L) \ for (quasi) 
one dimensional disordered samples can be calculated an- 
alytically JlTj . It indeed has a log-normal form for L ^> £, 
but all the positive moments (v n ) are dominated by rare 
events in such a way that (v n ) oc exp— cL/£, where the 
constant c is independent of the index n. 

To check if this property is shared by the level curva- 
tures we performed numerical simulations of the tight- 
binding Anderson model in one dimension. The curva- 
tures were calculated exactly from both eigenvalues and 
eigenvectors (see [?]), between 100 and 500 disorder re- 
alisations were used. We plot the numerical results for 
the logarithm of first, second and third moments of the 
curvature distribution versus L, see fig.l. The moments 
indeed decay roughly exponentially with system size, but 
the typical decay length decreases for higher moments 
in contrast to the behaviour typical for the quantity 

M-W(L)\. 

The second fact to be mentioned is that earlier numeri- 
cal investigations discovered a peculiar feature of the cur- 
vature distribution: the log-normal law is quite a good 
fit separately for the domain of curvatures larger and 
smaller than the most probable value (ln|if|) ~ 
but the parameters of the fit are slightly different for the 
two domains. At the same time, the statistics of the cor- 
relation function v = (— L)^f (L)\ is truly log-normal 
everywhere fL7| . 

These observations provide us with a motivation to 
consider the problem of the level curvature distribution in 
strongly localized regime on a more sound basis without 
invoking any additional assumptions. 

In the present paper we treat analytically the problem 
of the distribution of absolute value of the level curva- 
tures: 

V(K)=A(Y / 6(E-E n (0))5(K- \ K n |)\ (3) 

for a one-channel ring characterized by the Schrodinger 
equation: 



ring the mean free path is of the same order as the local- 
ization length. As we shall see, the negative moments of 
the curvature distribution Eq.(j|) can be found explicitly 
for rings with L » 1, i.e. in strong localization regime. 
This allows us to verify analytically the log-normal na- 
ture of the ensuing distribution at small curvatures. 

To address the problem of level curvatures in the most 
direct way we follow the method by Dorokhov [lq| and 
Kolokolov Jl9| who calculated the absolute value of the 
persistent current j n =| |. To this end, let us asso- 
ciate with any point x; —L < x < L within the sample 



a vector V(x) 



_ ( v+(x) 



with components: v±(x) = 



7J_ (x) 

±(d^/dx ± ik'fy) exp (^fikx) and consider a 2 x 2 trans- 
fer matrix T relating the value V(x) to "initial" value 
V(-L) in the following way: V(ar) = t(x,-L)V(-L). 
Due to the current conservation the transfer matrix can 
be parametrized as 



T(x,-L) = 



coshlV" sinhlV 3 
sinhTe"'' 3 coshre"'". 



(5) 



Here a(x, E),(3(x, E) and T(x, E) are real functions to be 
determined. The periodic boundary condition Eq.(|l|) can 
be written in terms of the transfer matrix T = T(L, —L) 
as det (f e 2( - lkLa ^ - e i( M = 0, where a z stands for the 

Pauli matrix. It is convenient to rewrite this condition 
in terms of the T- matrix elements as: 

f(E m , L) = cos 4>; f(E, x) = coshTcos (a + 2kx) (6) 

This equation determines the set of energy levels 
E m (<p). The following identities can be immediately in- 
ferred from this fact: 

HE ~ E n {4>)) = 8 (f(E, L) - cos <t>) \df/dE\ , (7) 



d 2 E n 



= 0) 



dE 



(E = E n ) 



(8) 



As a result one can write down the moments of the 
distribution Eq.(H) in the following form: 



M„ 



K n T(K) 



1 



A'- 



*(/-!) 



Of 



dE 



-71+1 \ 



(9) 



-— I + U(x)-E)y(x) = 



(4) 



with the boundary condition Eq. (jl|). Here U (x) is a white 
noise potential (U(x)U(x')) = D5(x — x') which is con- 
sidered to be weak: I = Ak 2 /D 3> fc _1 , with k being the 
Fermi momentum related to the energy E as k 2 = E and 
I standing for the mean free path. For the one-channel 



When performing the disorder averaging it is conve- 
nient to get rid of the 8— functions in Eq.(^) by aver- 
aging the corresponding expression over the ensemble 
of samples with slightly fluctuating sample lengths L ± 
8L; k~ x 8L <c min(i, L) |l9|]. When doing this we take 
into account that the functions T(x, E),a(x, E) change 
very slowly on the scale 8x ~ fc" 1 , i.e. (dT/dx, da/dx) <C 
t As a consequence there is typically only one solution 



2 



M n 



1 



1 1 






— tanhr 




\ | sinhT | 


dE 





(10) 



The disorder averaging in expressions of such a type 
can be performed by employing the functional integral 
method suggested by Kolokolov jli]]. For accomplishing 
such a calculation it is important to express the quan- 
tity to be averaged in terms of elements of the matrices 
T = T(L, — L) and dT/dE without involving complex 
conjugation. To find such a representation let us con- 
sider the auxiliary quantity: 



A = (0, 1) Ts's+T 



(11) 



where s are familiar lowering and raising operators for 
the spin 1/2. Taking the parametrization Eq.(^J) into ac- 
count we find that A — cosh 2 Te~ 2la and correspondingly 



dA/dE dT 



2A 



dE 



tanhT- 



. da dT 

'dE = dE 



tanhr 1 + O 



Lk 



(12) 



The right hand side of this expression is just the impor- 
tant part of the combination appearing in the expression 
Eq. ( |l0| ) under the sign of averaging. 

From the other hand the same quantity can be 

written in terms of the matrix dT/dE when differentiat- 
ing Eq.fjll]) with respect to E. 

The matrix dT/dE itself can be determined from the 
following exact relation between the solution ^>(x) of the 
initial Schrodinger equation and its derivative with re- 
spect to the energy $(a;) = d^(x)/dE. 



Cl 



dy 



CO 



(13) 



where c = ($ d^/dx - * d$/dx) | x =_£, and Ci = 
$(-£)/*(-£). 

We can write the relation: 



dT 
dE 



V(-£) = F(L) - TF(-L), 



(14) 



where F(x) 



dV(x 
dE 



The components f± of the vector F 



are represented as integral functionals of the fields v±(x), 
the latter fields having boundary values v±(— L) at the 
point x = —L. 

Combining all this facts together one can find the ap- 
propriate representation for the quantity dA J^ E ■ Skip- 
ping the cumbersome intermediate steps in favour of pre- 
senting the final expression we find finally: 



2A 



of the equation Eq.(^) in the interval Sx — ir/k. We dA/dE 
also found that the quantities Lk^ 1 are of the lower 
order 0(1/ (Lk)) when compared with ^ and can be 
safely neglected. This gives the possibility to rewrite the 
expression Eq.(^|) in the form: 



dy 



-l v+(y)v-(y) J_ 



v+(yi)v-(y 1 )dy 1 , (15) 



where the components v±(x) of the field V(x) are given 

by V(a;) = T(x,—L) f ^ . In the expression Eq.(|l5|) 

we kept only the leading terms with respect to 1/ Lk and 
omitted also all terms containing the fast oscillating fac- 
tors e ±2lkL . 

The following comment is appropriate here. Gener- 
ally speaking, the functions v±(x) are random and can 
take any complex values. However, they are not inde- 
pendent from each other and can be chosen to satisfy the 
constraint v+(x) = w* (x). When using the Kolokolov's 
method of averaging the random fields, v + (x) and w_(x) 
should be analytically continued on the surface deter- 
mined by this constraint to provide the convergence of 
the corresponding path integral. Apart from this fact, in 
our case one can use the same constraint to argue that 
da/dE is small in comparison with tanh T dT/dE by the 
factor l/Lk which we intensively used to simplify expres- 
sions give n above, e.g. Eq. (fl2"|) . 

Eq.(^o|) combined with Eqs.(|l2|- ^5|) provides a repre- 
sentation which is used as a starting point for employ- 
ing the Kolokolov's approach. Indeed, we expressed the 
moments M n in terms of the elements of the T-matrix 
without complex conjugation: the components of v±(x) 
were related to the T-matrix above and | sinh T \ can be 
represented in terms of the T-matrix in the same way as 
in p|. 

After the set of manipulations identical to those used in 
p9| we represent the moments in the form of the following 
path integral: 



-oo £(_L) =( j'; £(L)=<r 



V^e 2 x 



/ e-^dx J e^dy 

L -L 



-n+1 



(16) 



where a = 2/1. 

In general, only the first moment (n = 1) related to the 
average persistent current can be explicitly evaluated for 
arbitrary relation between the mean free path I and the 
ring circumference 2L Jlgjjul . 

However, in the limit 2L 3> I we can extract the leading 
contribution to the path integral Eq.(|l6|) for a large, but 
limited number of negative moments satisfying | n \ < 
aL/2. To see this let us use the notation m = —n so that 
the twofold integral in Eq.(|l6|) is raised to the positive 
power m + 1. It gives a possibility to treat the path 
integral in Eq.(|l6|) as the sum of 2(m + l)-fold integrals 
from the corresponding matrix elements. To illustrate 
the structure of the terms we write down explicitly the 
simplest one: 
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M = 
1 



^ e B « I-L dx I-L d y x 

L)H 



e -(L-x)H t 2 e -(x-y)HJ_ e -(y- 

t 2 



(17) 



where we introduced the new variable i = (2/a)e 

in terms of which the Hamiltonian operator TL and the 

scalar product (f\g) read K = | (t 2 <9 2 + i<9 t — i 2 ) and 

(/Iff) = Jo" f correspondingly. 

Higher moments will contain similar integrals over vari- 
ables — L < xi, yi, X2, J/2, ■■■■x-m, Um < L. Each integral 
of this type can be evaluated analytically if we expand 
powers of t in terms of the eigenfunctions of the opera- 
tor 7i which are the modified Bessel functions:7i A p (i) = 
—a/8p 2 K p (t). The functions with imaginary indices 
p = iv form a complete orthogonal set suitable for expan- 
sion. In this way we find that the leading contribution to 
the path integral Eq.(|l6|)m the limit L I corresponds 
to the configuration: yi w yi rj ... « y m = —L; x\ ~ 
... w x m = L when we can effectively write: that 



m+l 



Km + l)!] 2 a - 2 ( m + 1 ) e -( m + 1 )«( L ) e ( m + 1 )?(-- L ) 



(18) 



for < m <C aL/2. For larger m the contribution from 
omitted terms becomes comparable with that given by 
Eq.(|l8[) due to large combinatoric factors. 

The moments M_ m in this approximation is equal to: 



M_ m = 



1 



,-aL/A 



(19) 



When restoring from these moments the asymptotic be- 
haviour of the curvature distribution the factorial coeffi- 
cients can be omitted in view of the condition m -C aL/2. 
This results in the log-normal curvature distribution: 



V(hxK) ~ (47raL) 



(20) 



which is valid in the strongly localized limit aL 3> 1 in- 
side the parametricallly large domain of small curvatures 
-(aL) 2 < In A' < -aL. 

In conclusion, we demonstrated explicitly that the 
statistics of level curvatures in ID disordered systems is 
log-normal inside a parametrically large domain of cur- 
vatures smaller than the typical value. At the same time, 
interesting issues of explaining an unexpected difference 
in behaviour of the positive moments of the curvature 
and that of the eigenfunction correlator, see fi.g.1, as well 
as the mentioned asymmetry of the curvature distribu- 
tion remain open. 
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FIG. 1. 1st, 2nd and 3rd moment of the curvature dis- 
tribution as function of the system size as obtained from 
numerical simulations of the tight -binding Anderson model 
in one dimension. The disorder parameter w (width of the 
box-distribution of the diagonal matrix elements) is w — 2. 
Typical error bars are comparable to the remaining fluctua- 
tions on top of the exponential decay. 
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